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ABSTRACT 

We present a spectrophotometric study of short-term optical variability in the qui- 
escent black hole X-ray transient V404 Cyg. This includes two nights of high time- 
resolution Ha spectroscopy with which we resolve much of the time-variability, and 
a further six nights of archival spectroscopy with lower time-resolution but higher 
spectral-resolution. We find significant variability in most of the data considered, with 
both the Ha line and the continuum often varying in a correlated way. This includes 
both dramatic flares lasting a few hours in which the line flux nearly doubles and 
lower-level flickering. The strongest flares involve development of asymmetry in the 
line profile, with the red wing usually strongest independent of orbital phase. It is 
unclear why this is the case, but we discuss several possible explanations. We consider 
the energetics of the flares and compare with plausible models including chromospheric 
activity on the companion star, local magnetic reconnection events within the disc and 
a varying irradiation from close to the black hole. Based on the line profile changes 
during the flares, we conclude that the most likely origin for the variability is variable 
photoionisation by the central source, although local flares within the disc cannot be 
ruled out. 
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1 INTRODUCTION 

Quiescent black hole X-ray transients (BHXRTs) provide 
the best evidence we have for the existence of stellar mass 
black holes. These X-ray binaries contain a relatively low- 
mass star (usually, but not always, < IMq) accreting onto 
a likely black hole via Roche lobe overflow and an accre- 
tion disc. In their quiescent state the accretion flow becomes 
extremely faint and so the companion star can be directly 
observed. It is the companion, and in particular its radial 
velocity variations, which provide the key to measuring the 
black hole mass (Charles 1998). 

One of the most secure cases for a black hole is for 
V404 Cyg, which contains a compact object with a mass 
function of 6.08±0.06 M (Casares, Charles & Naylor 1992), 
twice the usually accepted upper limit to a neutron star mass 
(Rhoades & RufRni 1974). The best estimate for the actual 
mass is 12±2 Mq (Shahbaz et al. 1994) and this system likely 
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contains the most massive stellar mass black hole for which 
we have dynamical evidence. It has the longest orbital period 
of the BHXRTs, at 6.5 days, with a KOIV companion star 
(Casares et al. 1993). V404 Cyg is the most luminous of 
these systems in quiescence with Lx ~ 10 33 — 10 34 ergs -1 , 
~ 1000 x brighter than the prototypical BHXRT, A 0620- 
00 (Garcia et al. 2001). These extremes are unlikely to be 
coincidental and the high Lx may be a consequence of the 
long orbital period, the evolved donor star, and/or the high 
mass of the black hole. 

The high quiescent Lx indicates that even in its quies- 
cent state, V404 Cyg is not completely dormant, but that 
some accretion continues. This is also indicated by the sub- 
stantial short-term (sub-orbital) variability that V404 Cyg 
exhibits. In the optical this has an amplitude of 0.1-0.2 mag 
(Wagner et al. 1992), at least partly due to a 6hr quasi- 
periodic oscillation (Casares et al. 1993; Pavlenko et al. 
1996). Short timescale infrared variations also appear to be 
present (Sanwal et al. 1996) and the Ha line profiles change 
significantly on short timescales (Casares & Charles 1992; 
Casares et al. 1993). At X-ray energies, Wagner et al. (1994) 
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saw changes of a factor of 10 on timescales < 0.5 day with 
ROSAT and more recently Garcia et al. (2000) presented a 
Chandra lightcurve exhibiting a factor of ~ 2 variability in a 
few ksec. Finally V404 Cyg remains a variable radio source 
on timescales of days: Hjellming et al. (2000) reported typ- 
ical variability in flux of 0. 1-0.8 mjy with a maximum qui- 
escent flux of 1.5 mjy. 

V404 Cyg is not the only BHXRT to exhibit variability 
in quiescence. Haswell (1992) discovered significant optical 
flaring in A 0620-00 on timescales of minutes. Zurita et al. 
(in preparation) have found optical photometric variability 
in other systems: in addition to V404 Cyg and A 0620-00, 
they also identify strong variability in GRO J0422+32 and 
also in the neutron star system Cen X-4. This kind of be- 
haviour thus seems to be common in quiescent BHXRTs. 

It is important to understand these variations as the 
properties of the variability can yield important clues about 
the nature of the quiescent accretion flow and provide an 
observational test of models for this flow. For example, ad- 
vection dominated accretion flow (ADAF) models for V404 
Cyg have predicted that the non-stellar optical continuum 
is dominated by synchrotron emission from very close to the 
black hole (Narayan, Barret & McClintock 1997; Quataert & 
Narayan 1999). One would most naturally expect the ADAF 
temporal variability to be different to that of the emission 
lines which show a double peaked profile, indicative of an 
origin in the outer disc, although there is also a possibility 
that the lines could be photoionised by the ADAF. Alterna- 
tively, both line and non-stellar continuum emission could 
originate in the outer disc, with variability likely associated 
with local magnetic reconnection events. One would expect 
in this case rapid and kinematically localized enhancements 
in the line emission, although the flares might be broadened 
by motions within the flaring material. Another important 
source of variability may be the companion star. This is a 
K0 IV star, tidally locked in a 6.5 day orbit and likely similar 
to active stars in RS CVn binaries. 

To investigate these issues, in 1999 July we obtained 
optical spectroscopy at high time-resolution so as to reveal 
the spectral signature of the short-term variability. We ob- 
tained simultaneous photometry to ensure a reliable flux 
calibration; hence we can accurately study both line and 
continuum variations. The low spectral resolution also pro- 
vides us with kinematic resolution of emission line varia- 
tions. To place these observations in a broader context we 
include longer term photometric monitoring observations 
taken during 1999. We also re-analysed high spectral resolu- 
tion data from a spectrophotometric campaign on V404 Cyg 
performed in 1992 July. A preliminary analysis of the lat- 
ter data appeared in Zurita, Casares & Charles (2000). We 
describe the observations in Section ^| Section g presents 
the average spectrum of V404 Cyg from 1999 together with 
that of the nearby blended star. In Section l| we compare 
our line and continuum lightcurves. Section ^ examines the 
spectrum of the variability in more detail and in Section ^ 
we discuss our results in the light of current models for qui- 
escent BHXRTs. 



2 OBSERVATIONS 

We observed V404 Cyg from the Observatorio del Roque 
de los Muchachos on La Palma on 1992 July 5-6 and 8-11 
and 1999 July 6-8 using the 1.0 m Jacobus Kapteyn Tele- 
scope (JKT) for photometry and the 4.2 m William Herschel 
Telescope (WHT) for spectroscopy. Conditions on all nights 
were good, with little or no dust or cirrus and seeing typi- 
cally ~ 1 arcsec. Further photometry was obtained with the 
0.38 m telescope of the Crimean Astrophysical Observatory 
(CrAO) and the 1.25 m and 0.6 m telescopes of the Stern- 
berg Astronomical Institute (SAI) from 1999 September 11 
to November 19. A full log of observations used is given in 
Table 0. 



2.1 JKT Photometry 

To ensure accurate photometric calibration of the WHT 
spectroscopy, we obtained simultaneous R band photome- 
try. In 1999 we used the SITe2 CCD camera on the JKT. 
We used integration times of 60 s on both nights, obtaining 
a temporal resolution of ~80s. All images were corrected 
for bias and flat-fielded in the standard way using iraf. Al- 
though the seeing conditions were good, we were not able 
to cleanly separate the contribution of the target and its 
nearby contaminating star so we performed simple aper- 
ture photometry around both stars. The contribution of the 
blended star, assumed constant, was then subtracted assum- 
ing R — 17.52 (Casares et al. 1993). We used an aperture of 
9 pixel radius (~3 arcsec) for our object and several compar- 
ison stars within the field of view, which were checked for 
variability during each night and during the entire data set. 
We also explored any correlation between the target differ- 
ential variations and the seeing, but found none. The 1992 
data were obtained with the GEC CCD camera on the same 
telescope. The typical integration time was 200 s. Aperture 
photometry was performed with an aperture of 11 pixel ra- 
dius (~3. 6 arcsec). Further details of these observations can 
be found in Pavlenko et al. (1996). 



2.2 WHT Spectroscopy 

We obtained the 1999 time-resolved spectrophotometry us- 
ing the ISIS dual-arm spectrograph on the WHT. To max- 
imise efficiency and minimise readout time and noise, we 
used the single red-arm mode with the relatively low- 
resolution R316R grating and the TEK4 CCD. Cycle times 
were 203 s on 1999 July 6-7 and 263 s on 1999 July 7-8. A 
wide slit (4 arcsec) was used to maximise photometric accu- 
racy, so our spectral resolution was determined by the seeing 
(~ 1.0 arcsec), and was typically ~ 4 A. Bias correction and 
flat fielding were performed using standard iraf techniques. 

Observations of V404 Cyg are made difficult by the 
presence of a fainter star 1.5 arcsec to the north and usu- 
ally blended with it (Udalski & Kaluzny 1991; Casares et 
al. 1993). Because of this, we chose a slit alignment pass- 
ing through the line of centres of the two stars. Fortunately 
this also coincided with a third, isolated, relatively bright 
star. Hence we could use the isolated star to define the spa- 
tial profile along the slit, and then deblend the spectra of 
V404 Cyg and the nearby star, following the method sum- 
marised in Hynes et al. (1999) and described in more detail 
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Table 1. Log of spectroscopic and photometric observations of V404 Cyg from 1992 (La Palma) and 1999 (La Palma, CrAO and SAI). 



Date 


Telescope 


UT Range 


Exposures 


Number of 


Wavelength 


Resolution 








Time (s) 


Exposures 


range (A) 


(A) 


La Palma Spectroscopy 














1992 Jul 5 


WHT 


00:51-04:43 


1800 


8 


6200-6600 


0.5 


1992 Jul 6 


WHT 


00:38-05:04 


1800 


8 


6200-6600 


0.5 


1992 Jul 8 


WHT 


00:33-04:58 


1800 


10 


6200-6600 


0.5 


1992 Jul 9 


WHT 


00:56-04:24 


1800 


7 


6200-6600 


0.5 


1992 Jul 10 


WHT 


00:20-04:53 


1800 


9 


6200-6600 


0.5 


1992 Jul 11 


WHT 


23:46-03:43 


1800 


8 


6200-6600 


0.5 


1999 Jul 6 


WHT 


21:27-04:02 


180 


98 


6000-7400 


4 


1999 Jul 7 


WHT 


21:30-04:04 


240 


83 


6000-7400 


4 


La Palma Photometry 














1992 Jul 5 


JKT 


01:09-05:15 


200 


33 






1992 Jul 6 


JKT 


00:34-05:17 


200 


50 






1992 Jul 8 


JKT 


00:28-05:11 


200 


63 






1992 Jul 9 


JKT 


00:56-04:31 


200 


44 






1992 Jul 10 


JKT 


00:13-04:40 


200 


49 






1992 Jul 11 


JKT 


23:40-04:10 


200 


45 






1999 Jul 6 


JKT 


22:03-03:56 


60 


205 






1999 Jul 7 


JKT 


21:54-03:57 


60 


259 
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1999 Sep 17 


38 cm 


20:12-00:01 


180 


43 






1999 Sep 18 


38 cm 


19:26-20:24 


180 


11 






1999 Sep 21 


125 cm 


22:31-01:34 


90 


96 






1999 Sep 22 


125 cm 


22:43-00:45 


90 


56 






1999 Sep 23 


125 cm 


23:26-02:09 


90 


73 






1999 Oct 16 


38 cm 


20:52-22:04 


180 


23 






1999 Nov 05 


38 cm 


20:15-21:07 


200 


11 






1999 Nov 11 


125 cm 


19:21-19:26 


20 


2 






1999 Nov 14 


125 cm 


17:59-18:14 


90 


7 






1999 Nov 19 


60 cm 


18:28-18:33 


120 


3 







by Hynes (1999,2001). For these data, we used the isolated 
star to construct a wavelength independent semi-analytical 
spatial profile for each spectrum. This comprised a Voigt 
profile with a numerical correction determined for the pro- 
file wings. The profile deblending optimal extraction method 
does reject most cosmic rays and bad pixels, but a few re- 
mained; these were corrected by hand. 

Wavelength calibration was interpolated between con- 
temporaneous exposures of a copper-neon/copper-argon arc 
lamp. The dominant residual error was caused by the uncer- 
tain positioning of the star within the wide slit. The sky lines 
could not be used to correct the wavelength solution because 
these fill the slit and the stars positioning relative to these 
is also uncertain; hence we must use spectral features of the 
objects. Unfortunately the template star shows no strong 
spectral features in this region. We experimented with us- 
ing both the Ha absorption feature in the faint blended star 
and various stellar absorption features in V404 Cyg (over 
AA6380-6520 A) and found that the latter gave less scatter 
in the corrections, which were calculated by cross-correlating 
the individual spectra in this region with an average. Of 
course the photospheric features do move, so this correction 
represents a transformation into the companion rest frame. 
We then use the known radial velocity curve of Casares & 



Charles (1994) to transform into an inertial frame. Initial 
flux calibration was done relative to the template star. This 
was in turn calibrated relative to the spectrophotometric 
standard Cyg OB2, No. 9 (Massey et al. 1988). Finally we 
performed synthetic R band photometry of the sum of cal- 
ibrated spectra of V404 Cyg and the blended star and slit 
losses were corrected by comparing with the simultaneous 
JKT photometry. 

V404 Cyg was also observed in the Ha region with the 
ISIS spectrograph on the WHT on the nights 1992 July 5-6 
and 8-11. A total of 50 spectra were obtained using indi- 
vidual exposures times of 1800 s. The main aim of these 
observations was to measure the rotational broadening of 
the companion star, so a narrow (0.8arcsec) slit was used in 
combination with a 1200 line mm -1 grating to yield a high 
spectral resolution of 0.5 A (FWHM). The wavelength range 
covered was ~ AA6200-6600 A. The spectra were calibrated 
using the flux standard BD+17°4708 (Oke & Gunn 1983). 
Further details of these observations can be found in Casares 
& Charles (1994). As the slit was narrower, the slit loss cor- 
rection is more dramatic than in the 1999 observations. We 
made the correction, as is explained in the first part of this 
section, by comparing the simultaneous JKT photometric 
fluxes and the spectroscopic fluxes. 
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Figure 1. Deblended spectra of V404 Cyg and the blended star. 
All the spectra from 1999 July have been averaged. The deblend- 
ing algorithm has clearly separated the spectra cleanly. Note the 
contrast between He* emission in V404 Cyg and absorption in the 
blended star. 



2.3 Contemporaneous Photometry 

As part of an ongoing campaign to monitor V404 Cyg 
in quiescence, observations were carried out at the 0.38 m 
Cassegrain telescope of the Crimean Astrophysical Obser- 
vatory (CrAO) and the 1.25 m and 0.6 m telescopes of the 
Sternberg Astronomical Institute (SAI) in a broad R + I 
system, using SBIG ST-7 cameras. The images were dark- 
subtracted, flat-fielded and analysed with a profile/aperture 
photometry package developed by Goransky (SAI). 



3 AVERAGE SPECTRA 

In Fig. ^ we show the average spectra of V404 Cyg and the 
nearby blended star from 1999. To our knowledge, this is the 
first time that the spectra of the two stars have been sep- 
arated; most spectroscopic observations either yield a com- 
bined spectrum, or align the slit perpendicular to the line 
of centres to minimise the contribution from the blended 
star. Contamination by this star will distort the spectrum 
and must be considered (e.g. Shahbaz et al. 1996) so it is 
worthwhile to present the separated spectra for comparison. 

Photometry of the blended star yielded B — V = 1.69, 
V — R = 1.38 (Casares et al. 1993). For different trial redden- 
ings, the best agreement with colours of main sequence stars 
comes for a K4-5 star with E(B — V) ~ 0.5. The reddening 
line in a B - y vs. V — R plot is almost parallel to the main- 
sequence for spectral types earlier than K, however, so more 
reddened earlier types are also possible. The only strong fea- 
ture visible in the blended star's spectrum is Ha absorption. 
This is stronger than expected for a late K star and instead 
favours an F type identification; a higher resolution spec- 
trum would be needed for a more precise classification. A 
weak 6495 A absorption line also appears to be present; this 
feature can be present in F stars (see e.g. the templates used 
by Orosz & Bailyn 1997), so is consistent with this interpre- 
tation. For the photometric colours to agree with an F star 
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Figure 2. Phase-folded photometric lightcurves from CrAO and 
SAI (solid circles) and the JKT (open circles). All lightcurves have 
been binned in groups of three to improve signal to noise, and 
the R band JKT points have been shifted upwards by 0.38 mag 
to bring them into approximate agreement with the R + I band 
observations. 



requires E(B - V) ~ 1.3 ± 0.1, or A v ~ 4.0 ± 0.3. This im- 
plies a dereddened brightness of V ~ 14.9 and a distance of 
2.0 ± 0.5 kpc for a main sequence star. Both reddening and 
distance are similar to the values for V404 Cyg (Casares 
et al. 1993) so this is a plausible solution. Further support 
for a similar reddening to V404 Cyg is the presence of the 
6284 A diffuse interstellar band with a similar strength in 
both objects. If the blended star is an F star then we ex- 
pect it to have no significant CO features in the infrared. 
Consequently the infrared disc contribution will be reduced 
further below the 14 per cent upper limit of Shahbaz et al. 
(1996). 

For the following analysis we will require the Ha 
line luminosity. We observe an average line flux of 4.8 x 
10~ 15 erg cm -2 s _1 , and Casares et al. (1993) estimate an 
optical extinction of Ay = 4.0. Assuming the gas to dust 
scaling of Bohlin, Savage & Drake (1978), this extinction cor- 
responds to a neutral hydrogen column of ~ 7 x 10 21 cm -2 , 
which is also consistent with that inferred from the Chandra 
spectra (Garcia et al. 2000). With Av = 4.0 and an average 
Galactic extinction curve (Seaton 1979) we infer an inter- 
stellar attenuation at Ha of a factor of ~ 20. A distance of 
3.5 kpc was adopted by Garcia et al. (2001) and to facili- 
tate comparison with their X-ray results we adopt the same 
value. Finally, to convert observed fluxes into luminosities we 
also assume isotropic emission. This is probably not correct 
for emission from a disc, but as V404 Cyg appears to be an 
intermediate inclination system (Shahbaz et al. 1994), the 
difference in the ratio of observed flux to luminosity between 
isotropic and disc emission is likely to be small, and this is 
acceptable for a crude approximation. We finally derive an 
average, dereddened Ha luminosity of 1.4 x 10 32 ergs _1 . 



4 LIGHTCURVES 
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Figure 3. Continuum light curves of V404 Cyg and the nearby blended star, and Ha emission light curve of V404 Cyg from the 1999 
spectrophotometry on July 6—7 (left) and 7—8 (right). Continuum measurements are based on a weighted average across the spectral 
range covered, with lines masked out and with higher weight given where the signal-to-noise ratio is highest. Note that the vertical range 
of the Ha plot is much larger than for the continuum light curves. The blended star is fainter than V404 Cyg, so provides an upper limit 
to the statistical and systematic errors expected in the lightcurve of Ha. The notation for dividing up the lightcurve is Q: 'quiescent 
state', F: isolated flickering regions, somewhat brighter than Q and fR, fP and fD denote the flare rise, peak and decline respectively. 



4.1 The orbital lightcurve 

In Fig. ^ we present phase-folded photometric lightcurves 
from CrAO, SAI and the JKT in 1999. Scales have been 
chosen to facilitate comparison with Fig. 1 of Pavlenko et al. 
(1996). As with the 1992 data described there, the lightcurve 
is dominated by the double-humped ellipsoidal modulation, 
but superposed upon it are shorter term variations visible 
within a single night. The 1999 data were obtained in a red- 
der bandpass (R + I) than those obtained in 1992, so are 
not directly comparable, but it does appear that the am- 
plitude of both ellipsoidal and short-term variations is less. 
Apparent changes in the ellipsoidal lightcurve of V404 Cyg 
based on the monitoring campaign have already been noted 
by Pavlenko et al. (2001). 

In the earlier data the short term variations appeared 
quasi-periodic (Pavlenko et al. 1996). We performed a pe- 
riod search on the 1999 data to test if this was still the case. 
Before doing the analysis we removed the 6.4714 d orbital 
modulation and its first two harmonics. The periodogram 
of the residuals showed several peaks and it is difficult to 
choose between them, or even prove that the variations are 



quasi-periodic. Formally the peak with the highest signifi- 
cance has period (0.25404 ± 0.00005) d, ~ 6.1 hr with mean 
amplitude ~ 0.1 mag. This is close to the period found in 
the earlier data, but other periods are possible and the case 
is inconclusive. 



4.2 The high temporal resolution spectra 

Light curves from the 1999 high temporal resolution spectra 
are shown in Fig. ^| On both nights we find the variability 
is dominated by a single large flare, although the profiles of 
these flares are complex. On the second night the flare is 
clearly seen in the continuum as well, with the same struc- 
ture. The amplitude of the line flare is much larger than that 
in the continuum, but this is likely because the continuum 
is heavily diluted by light from the mass donor star; Casares 
et al. (1993) estimate 90 percent of the light in the R band is 
from the mass donor. After correction for this, the fractional 
amplitudes of the flares in line and inferred accretion flow 
continuum are similar. On the first night, the continuum 
does not reflect the large flare, although finer structure does 
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appear common to the lightcurves and a weak continuum 
flare does coincide with the peak of the Ha flare. 

To test if the flaring is simultaneous, or if one lags 
the other, we calculated cross-correlation functions between 
them. The lightcurves are reasonably well sampled so we use 
the interpolation cross correlation function (Gaskell & Pe- 
terson 1987; White & Peterson 1994). We also calculated dis- 
crete correlation functions (Edelson & Krolik 1988) and the 
results were almost identical. Our CCFs for the two nights 
are shown in Fig. ^. As noted above, the correlation between 
line and continuum does not show up well on the first night, 
and at 3 cr confidence is not significant. On the second night a 
clear peak is seen, significant at the 3 a level. The peak shows 
no large lag between line and continuum, but does appear 
asymmetric. This could correspond to line variations lasting 
longer than those in the continuum, perhaps due to a sig- 
nificant recombination timescale. The lack of correlation on 
the first night, and the strange structure of the CCF on the 
second night are at least partly due to low-frequency vari- 
ations, in particular the absence of large continuum flares 
when large Ha flares occur. On the first night there is no 
prominent continuum flare, but some shorter timescale vari- 
ations do appear correlated. On the second night only the 
peaks of the large flare are well reproduced by the con- 
tinuum, not the broader base. To further explore this, we 
removed the poorly correlated low-frequency variations by 
applying a high-pass Fourier filter to the lightcurves with a 
cut-off at lhr -1 . The CCFs of the filtered lightcurves are 
shown in the lower panel of Fig. ^. On the first night, there 
is still no prominent peak in the CCF, in spite of appar- 
ent similarities in the lightcurves. For the second night, the 
structure at negative lags has been removed, but the peak 
remains prominent and definitely appears significant. We es- 
timate a delay of the line with respect to the continuum from 
the filtered lightcurves of (70 ±40) s, consistent with zero at 
the 2 a level. The error estimate follows Gaskell & Peterson 
(1987) and is defined for uniformly sampled data. Our data 
are close to, but not perfectly, uniform. Hence we believe this 
error estimate is not too far wrong. We also calculated dis- 
crete CCFs (Edelson & Krolik 1988) to check that the peak 
was not being pulled towards zero lag by correlated noise 
in the line and continuum lightcurves; this method elim- 
inates correlated noise by only including the contribution 
of line and continuum from different spectra. The discrete 
CCF actually yielded a slightly shorter lag of ~ 55 s. Since 
the measured delay is not significant at the 3 a level, and 
is below the time resolution used, it is probably not real. 
The 3 a upper limit on the lag is then about 200 s. We note 
that important short timescales such as the light travel time 
across the binary (~ 30 s) are not resolved, so a smaller lag 
may be present; higher time-resolution observations would 
be required to test this, but would then require an 8-10 m 
class telescope. 

For the analysis of the line profiles in Section we di- 
vided the lightcurves up into segments, a convenient but ar- 
tificial division of continuously varying behaviour. We iden- 
tify the lowest region of each night with the 'quiescent state', 
Q. Isolated flickering regions, somewhat brighter than Q are 
labelled F. Within the flares we separate flare rise (fR), peak 
(fP) and decline (fD) segments. This nomenclature is sum- 
marised in Table 0. 

For the two major flares we estimated the peak flare 
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Figure 4. Cross-correlation functions between line and contin- 
uum from 1999 data. A positive lag would correspond to line 
variations lagging behind those in the continuum. Dotted lines in- 
dicate 3 a limits on expected coincidental correlations (almost the 
same on the two nights). Upper panel shows the cross-correlations 
of the raw light curves; lower panel is after applying a high-pass 
Fourier filter with a cut-off at 1 hr -1 . 



luminosity emitted in Ha and the total energy released in 
this line. We assumed that the lowest level reached on each 
night represents the base, persistent flux and integrate over 
the range actually observed. Since both flares extend beyond 
our coverage this is a lower limit to the total energy. For the 
first night we estimated a peak flare luminosity (at 3 kpc) of 
1.4 x 10 32 ergs -1 and a total observed Ha energy of ^ 9 x 
10 35 erg. For the second night the values are 1.5 x 10 32 ergs -1 
and <; 9 x 10 35 erg respectively. The peak luminosities of the 
flares are strikingly similar, although the total energies are 
likely different, since most of the decline of the flare on the 
first night was missed. 

4.3 The high spectral resolution spectra 

Light curves from the 1992 high spectral resolution spectra 
are shown in Fig. ^[ Here we can see that variability is also 
dominated by single flares with changes in flux up to a factor 
of 2 on time scales of 1-2 hr, although due to the poorer time 
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resolution, we cannot resolve any short-term structure. The 
largest flares are seen on the nights of 6, 9 and 11 July and 
both the amplitude of the events and the duration is similar 
to that found in 1999, seven years later. Generally the con- 
tinuum lightcurve is correlated with the Ha one, although 
the amplitude of continuum variations is ^ 50 percent of the 
Ha variability. 



5 SPECTRAL VARIABILITY 

Lightcurves are a good way to contrast line and continuum 
variability but information about velocities is lost. In this 
section we explore several methods to analyse the changes 
in the line morphology and constrain the binary region re- 
sponsible for the Ha variability. 



5.1 Trailed Spectra 

Our trailed spectra from 1999 are presented in Fig. []. It 
can be seen that the line profile varies with the flares, with 
strong red peaks seen at flare maxima. Other than this en- 
hancement of the red wing relative to the blue, we find no 
evidence of velocity changes, or new velocity components 
or lines appearing during the flares. There is an apparent 
broadening of the profiles in the trailed spectra during the 
flares, but this is not real, as can be seen by the profiles in 
Fig. \A it is due to the overall brightening of the profile. 



5.2 Line profiles from the high temporal 
resolution spectra 

We have constructed average line profiles from each of the 
lightcurve regions indicated in Fig. j3|. These are shown in the 
left-hand panel of Fig. Q During faint (Q) periods the Ha 
line is largely symmetric and exhibits a 'classical' double- 
peaked disc profile. During flickering (F) periods a similar 
symmetric disc profile is seen, but it is stronger. Profiles 
on the rise to a flare (fR) are mostly similar, the exception 
being the second fR period on the first night (the highest 
dashed line in Fig. (?], just before the peak of the flare.) In 
this case, a mild asymmetry had developed and the profile 
was clearly developing towards the fP form (see below). Ex- 
cluding this case, F and fR profiles are indistinguishable. To 
test for more subtle differences between Q and F/fR profiles 
we combined all of the profiles of each type (excluding the 
anomalous fR from the first night) and normalised the aver- 
ages to the same total line flux. These normalised profiles are 
plotted in the centre panel of Fig. |?]. On both nights there 
is a tendency for F/fR profiles to be broader than Q ones. 
This suggests that the brightening mainly involves higher 
velocity material. In principal an average profile could also 
be broadened by err ors in wavelength calibration which, as 
described in Section 2.2, is not very reliable. The rms scat- 
ter in the wavelength corrections was 0.5 A and 1.2 A on the 
first and second nights respectively corresponding to 20 and 
50kms _1 at Ha. To explain the observed broadening re- 
quires calibration errors of order 100-150 km s , however, 
so the broadening is likely real. The spectra at and after the 
peak of the flares (fP and fD) are qualitatively very differ- 
ent from those at lower levels. On both nights, a strong red 



Table 2. Summary of the 'states' used to divide up the 1999 data. 
These are identified by lightcurve behaviour, but define different 
line profiles as well. 



Abbrev. 


Type 


Line profile 


Q 


Quiescent 


Weak, symmetric 


F 


Flickering 


Medium, symmetric 


fR 


Flare rise 


Medium, usually symmetric 


fP 


Flare peak 


Strong, enhanced red wing 


fD 


Flare decline 


Strong, enhanced red wing 



peak developed during the flare; the asymmetry is most pro- 
nounced on the second night. The blue peak is clearly also 
enhanced, as can be seen in the difference profile (fP/fD - 
Q) shown in Fig. [?]. There are changes in the asymmetry 
during the flares, but these are different on each night and 
there is no repeatable difference between fP and fD profiles. 
The velocity of the red peak is ~ 200 - 300 kms" 1 . The 
characteristics of profiles in different phases are summarised 
in Table §. 

We also constructed continuum spectra for each of these 
periods to test if the spectrum becomes redder or bluer dur- 
ing the flares. We can find no convincing change in con- 
tinuum colour between Q and fP/fD periods, although the 
quality of difference or ratio spectra between these periods 
is very poor. 



5.3 The high spectral resolution spectra 

For the 1992 data, the time-resolution is poorer making it 
hard to define regions of the lightcurves in the same way 
as was done in Section These data do, however, bene- 
fit from a higher signal-to-noise ratio and better wavelength 
calibration. They are thus well suited to calculation of the 
root-mean-sqare (rms) spectrum, as is commonly employed 
in AGN studies (see e.g. Wandel, Peterson & Malkan 1999 
for definition). We find that the line core, especially the red 
core (regions confined within ±400 km s -1 ), is the main con- 
tributor to the line variability (see Fig. ^) on July 5, 10 and 
11, peaking at ~ +200 km s -1 , comparable to the peak ve- 
locities in 1999. On July 6, however, it seems to be the blue 
core that varies strongest. The other two nights do not show 
pronounced changes in the symmetry of the lines. We cannot 
find any narrow components associated with the expected 
velocity of the secondary star or the hot-spot/gas-stream. 
The line appears to be more symmetric and narrower in the 
low states, as found in the higher time-resolution spectra 
discussed above. 



6 DISCUSSION 

It is clear that significant variability is present in both the 
continuum and the Ha line. We note that with an average 
equivalent width of ~ 20 A, Ha will make a small contribu- 
tion to R band photometry 2 percent, so R band photo- 
metric variability should be dominated by the continuum. 
Possibly there is more than one kind of variability present, 
with different physical origins. Very prominent are the large 
flares that dominate each night's light curve. These last a 
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Figure 5. Continuum and Ha emission lightcurves of V404 Cyg in 1992 July 5—6 and 8—11. Continuum lightcurves are based on integrated 
(unweighted) flux within the bandpass. Note that the vertical range of the Ha plot is much larger than for the continuum light curves 
and that the scale is the same as in Fig. N. 
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Figure 6. Trailed spectra in the vicinity of Ha from the 1999 observations on July 6-7 (left) and 7-8 (right). White gaps indicate when 
calibration frames were taken. 
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Figure 7. Changes in line profiles for the 1999 observations. Left-hand panel: Calibrated fluxes for all profiles. Center panel: F/fR and 
Q profiles normalised to the same peak flux with all profiles of the same type from each night have been averaged together. Right-hand 
panel: Difference between average fP/fD and Q profiles for each night. 



few hours each and occur once per night, making it likely 
that these can be identified with individual events of the 
6hr QPO described by Pavlenko et al. (1996). In addition 
to these large extended flares, shorter timescale variability 
is present with a lower amplitude. While this may originate 
from the same mechanism responsible for the larger flares, 
it could also represent a different process. The timescales 
involved in this higher frequency variability are comparable 
to those found in shorter period BHXRTs (Haswell 1992; 
Zurita et al. in preparation). This minute timescale flick- 
ering may be typical of quiescent BHXRTs, but the more 
dramatic flares could be unique to V404 Cyg and associated 
with its longer orbital period, higher black hole mass, higher 
X-ray luminosity and/or some other factor. Alternatively, if 
the variability timescale scales with orbital period, then the 
variations seen in other BHXRTs correspond to the large 
flares in V404 Cyg, and the counterparts of the short-term 
variation seen in the latter have yet to be found. 

Several variability sites are possible. An obvious can- 
didate would be local magnetic flares (reconnection events) 
within or above the disc. Other possibilities include instabil- 
ities in the accretion stream, resulting in variability at the 
stream-impact point, chromospheric activity on the compan- 
ion star and instabilities in the inner accretion flow resulting 
in variable photoionisation by a central X-ray source. 

6.1 Kinematic evidence for the flare location 

An obvious diagnostic of the flare location is the line veloc- 
ities, as these should show a radial velocity modulation cor- 
responding to the flare site. The asymmetry seen in the flare 
profiles obv iously complicates matters; this is discussed in 
Section 6.5. We would expect this asymmetry to be applied 



relative to the motion of the flaring component, so the un- 
derlying radial velocity modulation should still be apparent. 



Fig. shows the peak velocities of observed flares together 
with several predicted radial velocity curves. The compan- 
ion star's curve is taken from Casares & Charles (1994). The 
compact object's curve is also derived from this assuming 
q — 0.06 (Casares & Charles 1994). For the stream impact 
point we calculate ballistic stream velocities corresponding 
to impacts at 0.5, 0.6, 0.7, 0.8 and 0.9 times the effective lobe 
radius (Eggleton 1983). We assume q = 0.06, Afj = 12 M Q 
and i = 56° (Shahbaz et al. 1994). These parameters obvi- 
ously have some uncertainty, but at least illustrate the range 
of velocities which should be present. The relative constancy 
of the flare velocities at the range of orbital phases observed 
is evidence against the flares originating in the companion 
star's corona, and a hotspot origin is similarly unlikely. In- 
stead the flares appear to originate in or above the disc. In 
fact the lack of scatter in the flare velocities is also difficult 
to reconcile with localised flares in the disc; these should be 
scattered across its surface, with a spread in peak veloci- 
ties comparable to the width of the double peaked profile. 
This uniformity of flare velocity, together with the apparent 
participation of the whole double peaked profile in the flare 
(clearly visible in the difference profiles in Fig. argues that 
the Ha flare emission is spread over the whole disc. If this is 
true then the flares most likely originate in photoionisation 
by variable X-rays. Unfortunately, given the limited number 
of flares observed, none of these arguments are statistically 
conclusive, so we now examine other evidence to support or 
reject this interpretation. 

6.2 An active companion star? 

During quiescence the optical emission of BHXRTs is domi- 
nated by light from the secondary; hence the observed vari- 
ability could be due to chromospheric flare activity on the 
companion star. High levels of activity can be best main- 
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Figure 8. RMS spectra for the six nights in the 1992 dataset. 
Vertical dashed lines at ±400 km s _1 mark velocities from the red 
and blue cores and wings. 



tained via tidal locking in tight binaries, as occurs in RS 
CVn systems, so we can expect the same behaviour in 
LMXBs. Comparisons with RS CVn stars have already been 
made by Bildsten & Rutledge (2000) to explain the qui- 
escent X-rays. They conclude that the X-ray luminosities 
seen in most quiescent BHXRTs are consistent with an ori- 
gin in coronal activity on the companion star (Lx up to 
10 32 ergs _1 is seen in RS CVn systems; Dempsey et al. 
1993), but explicitly exclude V404 Cyg from this as its X-ray 
luminosity is too high. Lasota (2000), however, goes further 
and argues that coronal X-ray emission is expected to be too 
faint to wholly explain quiescent X-rays in any BHXRTs. 
Chandra spectroscopy supports this, indicating that the X- 
ray spectra of BHXRTs in general, and V404 Cyg in partic- 
ular, are harder than those of stellar coronae (Kong et al., in 
preparation). It remains likely that the companion stars are 
active and that they do make some contribution to the X-ray 
emission, but this probably does not dominate. We should 
then also ask if activity in the chromosphere of the com- 
panion star could be responsible for the optical variability 




-600 I 

0.0 



1.0 

Orbital phase 



Figure 9. Velocities of flare peaks from 1992 (open circles) and 
1999 (filled circles). The solid line shows the radial velocity curve 
of the companion star, the dashed line is the curve for the com- 
pact object and the dotted lines are a selection of ballistic stream 
velocities at stream impact points corresponding to disc radii of 
0.5—0.9 times the lobe radius. See text for details. 



observed. Although it has been generally assumed that opti- 
cal flares on RS CVn stars are a rare phenomenon, very high 
levels of activity have been reported in a few of these sys- 
tems. For example, UX Ari is a double-lined spectroscopic 
binary with spectral type G5V/K0IV and 6.44 day orbital 
period. It exhibits chromospheric emission originating from 
the cooler component. Ha flares with peak line luminosities 
ranging between ~0.7 and ~ 1.7 x 10 s0 erg s -1 have been 
observed (Montes et al. 1996). These are of the same order 
of magnitude as the maximum Ha luminosity observed in 
the Ha flares of other RS CVn systems such as V711 Tau 
(Foing et al. 1994) and HK Lac (Catalano & Frasca 1994). 
The dereddened peak luminosities in the V404 Cyg flares 
are ~ 1.5 x 10 32 ergs~ 1 on 1999 July 6-8, two orders of 
magnitude higher than in these RS CVn systems and so ap- 
parently not consistent with an origin in the chromosphere 
of the companion star, if it is comparable to an RS CVn star. 
This is open to question, however, as the companion star to 
V404 Cyg may be a stripped giant (King 1993). In this case 
the internal structure will be very different to a 'normal' RS 
CVn star with the same orbital period and spectral type. Of 
course, given the line profile of the flares (Fig. ^), with an 
enhancement of both wings, the kinematics of the flares are 
also not consistent with an origin on the companion star, 
and overall we feel this interpretation can be rejected. 



6.3 Variable photoionisation from the central 
source? 

The most likely site for the Ha flares is then the disc. We 
next ask how are they powered? The most likely possibilities 
are i) variations in the central X-ray source producing a 
variable photoionisation, or ii) local flares within the disc or 
disc corona powering the flares directly. 

Support for the variable photoionisation interpretation 
is provided by Casares et al. (1993) who find Ha/H/3/H7 ra- 
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tios consistent with case B recombination and conclude that 
the emission lines are photoionised. To further test this, we 
estimate the photon flux in Ha and compare with Chandra 
observations reported by Garcia et al. (2001). The Ha lu- 
minosity estimated in Section tfl is 1.4 x 10 32 ergs -1 , corre- 
sponding to 4.6 x 10 43 phots . 

While a straightforward comparison with the observed 
Lx is possible it would not be very meaningful. Any pho- 
tons above 13.6 eV could ionise hydrogen, and the absorp- 
tion cross-section is highest for the lowest-energy photons; 
these will cause ionisations at lower optical depths where 
the recombination spectrum is less thermalised. Further, 
the X-ray spectrum has photon-index F ~ 2 (Garcia et al. 

2000) , so low-energy photons can be expected to be much 
more numerous than their more energetic counterparts. Con- 
sequently, if photoionisation is important in powering Ha 
emission at all, it will be dominated by the unobservable 
EUV photons, and to estimate these we must either adopt 
a model or extrapolate from the X-ray data. Extrapolating 
from the Chandra spectrum (r ~ 2, Lx = 4.9x 10 33 ergs" 1 ), 
we expect a luminosity of ~ 10 34 ergs _1 above 13.6 eV or 
8 x 10 43 phots -1 . This is extremely uncertain; for exam- 
ple modest changes in the photon-index of ±0.2 can change 
these numbers by a factor of two. We feel that at present 
there is no definitive model for the high-energy spectra of 
quiescent BHXRTs that is more reliable in the EUV than 
an extrapolation. Simple ADAF models do make definite 
spectral predictions for V404 Cyg; for example those pre- 
sented by Narayan et al. (1997) predict an EUV flux above 
a pure extrapolation of the X-ray power-law due to an EUV 
bump from inverse Compton scattering of synchrotron pho- 
tons. Theoretical doubts about the possible presence of out- 
flows and convection call these models into question, how- 
ever (Narayan & Yi 1994; Igumenshchev, Chen & Abramow- 
icz 1996; Blandford & Begelman 1999). Models of ADAFs 
with winds for V404 Cyg have been considered by Quataert 
& Narayan (1999). By suppressing the central accretion 
rate, these increase the importance of bremsstrahlung ra- 
diation from the outer flow relative to inverse Compton 
emission and consequently harden the spectrum. This gives 
a poorer fit to observations than a wind-free ADAF, al- 
though this can be rectified by changing assumptions about 
the microphysics. In convection dominated accretion flow 
(CDAF) models emission is also expected to be dominated 
by bremsstrahlung radiation from large radii in the flow and 
spectral predictions are rather similar to those of ADAFs 
with winds (Ball, Narayan & Quataert 2001). While the un- 
broken power-law of the original ADAF models is in good 
agreement with the spectral data, these latter models require 
some fine tuning of parameters to achieve a satisfactory fit. 
Completely different scenarios exist, for example X-ray emis- 
sion from a corona above a cold disc (Nayakshin & Svensson 

2001) but these make no spectral predictions at all. Finally, 
we should remember that the X-ray flux is known to vary 
by up to a factor of ten on timescales longer than the 10 ks 
Chandra observation (Wagner et al. 1994). Consequently, 
the numbers above are only order of magnitude estimates at 
best. 

Adopting this (with large uncertainties) we require that 
the disc trap of order 50 percent of the emitted EUV pho- 
tons, if one Ha photon is produced per photoionisation. This 
is an extremely simplistic assumption, however, as many of 



the photoionisations considered will produce rather ener- 
getic photoelectrons. These could go on to collisionally ionise 
another atom, or produce bremstrahlung photons that pro- 
duce more photoionisations. A more secure limit is provided 
by comparing the input and output energies; we expect that 
the reprocessed emission from a disc element must be less 
than the input ionising flux, and so: 

Lria = /1/2-kx 

where /1 is the fraction of the high-energy emission Lx 
which is input into the disc and f 2 is the fraction of the 
input energy emitted in Ha. /1 depends on both the ge- 
ometry and the local conditions, as any photons which are 
Compton scattered off the disc will escape without causing 
a photoionisation. We can estimate an upper limit on f 2 in 
the case B recombination limit (supported by observations 
of V404 Cyg; Casares et al. 1993). In this case all Lyman 
line and continuum emission is assumed to be reabsorbed 
and the Ha emission is 20-30 percent of the total Balmer 
and Paschen line and continuum emission for temperatures 
of 5000-20000 K (Osterbrock 1987). Hence f 2 ;S 0.3. In 
reality f 2 could be significantly below this limit, as ener- 
getic photoelectrons may lose energy by other processes be- 
sides collisional ionisation: scattering off other electrons and 
bremsstrahlung may both be significant and will tend to 
thermalise the input energy. This problem will be more se- 
vere for the higher energy photons which will tend to deposit 
energy at higher optical depths within the disc and ulti- 
mately give rise to thermal continuum emission. For f 2 S3 0.3 
we then require /1 J; 0.05, although due to the uncertainty 
in the incident ionising flux this number is rather uncertain. 

In principal, photoionisation could occur on the com- 
panion star or the disc. For V404 Cyg, however, the mass 
ratio is small (g = 0.06; Casares & Charles 1994) and the 
companion star will absorb S; 1 percent of the emitted flux 
from near the central source. For the disc, the calculation is 
straightforward for the simplest ADAF models with emis- 
sion from close to the black hole. For an isotropic point 
source illuminating a thin, but somewhat flared disc we 
expect a fraction ~ h/r to be absorbed; hence we require 
h/r <; 0.05 for an ADAF radiating from the centre. ADAFs 
with winds or CDAFs, however, are expected to emit most 
of their radiation from larger radii and hence can more ef- 
ficiently irradiate the disc. The difference is probably not 
huge, however, perhaps a factor of a few higher. Unfortu- 
nately, given the large uncertainties in extrapolating the X- 
ray spectrum into the EUV region and the variability in the 
X-ray flux, none of these possibilities can be ruled out and it 
appears that it is energetically possible for all of the advec- 
tive models considered to produce sufficient photoionisation 
on the disc to account for the observed Ha flares, provided 
that not too much of the input energy is thermalised at large 
optical depths. Equally if X-ray emission originates from a 
corona above the disc (Nayakshin & Svensson 2001) then a 
large fraction (~ 50 percent) will be incident upon the sur- 
face, so this model is also plausible, even with emission of a 
significant fraction of the energy in the continuum. 

6.4 Local magnetic reconnection events? 

An alternative interpretation of the flares would be that they 
originate in local magnetic reconnection events in or above 
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the disc. It is widely believed that a dynamo mechanism 
operates in accretion discs, driven by the strong shear pro- 
duced by differential rotation (see e.g. Tout & Pringle 1992) . 
There is limited observational evidence for this, but it is 
known that Balmer line emissivity in quiescent dwarf novae 
seems to follow the same power law with rotation rate (oc 
as in chromospherically active stars (Home 
& Saar 1991). Consequently it has been proposed that the 
mechanism powering the Ha emission lines is the same in 
both active stars and discs, a dynamo. The effect of the 
strong shear is that regions of oppositely directed magnetic 
fields can develop within the disc (Tajima & Gilden 1987). 
These fields can reconnect explosively (Haswell, Tajima & 
Sakai 1992) giving rise to a flare. We also expect flux tubes 
to rise from the accretion disc surface into its corona because 
of the magnetic buoyancy (Parker) instability as occurs in 
chromospherically active stars, so reconnection and flaring 
could originate there. 

In solar flares, it is commonly accepted that every Ha 
flare is accompanied by a soft X-ray event, hence the effects 
of X-ray irradiation have to be taken into account and the 
ionisation and excitation of hydrogen together with the ther- 
malisation of the background electrons will be two important 
processes. Simultaneous Ha/soft X-ray (0.05-2 keV) obser- 
vation s of solar flares sho w a Ha to soft X-ray energy ratio of 
~5% (Acton et al. 1982). The same ratio has been ob served 
in the chro mospherically active dMe star G1644AB (Doyle 
et al. 1987), where the total Ha flare energy is 6 percent 



of the X-ray energy in the the 0.05-2 keV range. Although 
in V404 Cyg no simultaneous observations are available, a 
factor of ~2 variability in the X-ray flux has been observed 



(Garcia et al. 200C), so the peak flare luminosity is of the 



same order of magnitude as the average luminosity, assum- 
ing similar spectra. We can estimate the total energy in 
the 0.05-2 keV range by extrapolating the observed X-ray 
power- law (r ~ 2) to be ~ 6 x 10 33 ergs"\ The observed 
luminosity in Ha is ~ 1.5 x 10 32 ergs _1 , that is ~2. 5 percent 
of the X-ray energy. This is a similar order of magnitude 
to that seen in solar and stellar flares, although obviously 
very uncertain since a large extrapolation of the X-ray spec- 
trum is involved. The X-ray observations could therefore be 
consistent with a magnetic reconnection scenario, assuming 
that the ratio of Ha to X-ray emission does not radically 
rescale for flares in or above a disc. Given the uncertainties, 
an additional X-ray component from close to the black hole 
is certainly still possible even with this interpretation of the 
variability. 

A major problem for this model is that the timescales 
involved, with flare rise times of a few hours, are rather 
short. The line profile enhancements and the large amplitude 
of the flares suggest that these are not localised events but 
that the whole disc brightens on this timescale. Given the 
long period (6.5 d) of V404 Cyg it is hard to see how this 
could happen, as viscous timescales in the disc will be much 
longer than this, and hence the whole disc should not be 
coupled on this timescale. 



6.5 Line profile variations 

The development of asymmetry in the line profile is a fur- 
ther interesting issue, especially as it is usually the red wing 
that is affected; only one night, 1992 July 6 shows preferen- 



tial enhancement of the blue wing. In addition to the peaks 
summarised in Fig. |^ the single night of spectroscopy from 
1990 July 20 presented by Casares & Charles (1992) shows 
a red flare at orbital phase ~ 0.0 and in 1991, strong red 
wings were seen on August 8 and 9, corresponding to phases 
~ 0.05 and ~ 0.20 (Casares et al. 1993). This suggests that 
the enhancement in the red wing is not dependent on orbital 
phase, so does not indicate kinematics of localised flares. It 
is also notable that it only occurs when the line is strongest; 
during fP, fD and high fR periods in the nomenclature of 
Section ^[ During F and low fR periods, the line profile is 
enhanced, but symmetrically. 

It has long been known that solar flares commonly ex- 
hibit a red asymmetry in the Ha profile (Ellison 1943; see 
Svestka 1962 for a more comprehensive dataset). The asym- 
metry abruptly increases at the onset of the flare, attains 
a maximum velocity of 40-100 km s -1 and decreases before 
the Ha line reaches its maximum (Ichimoto & Kurokawa 
1984). Early papers interpreted it as due to absorption of the 
blue wing of the emission profile, but Ichimoto & Kurokawa 
(1984) associate it with a downward motion in the flare chro- 
mospheric region. The analogue for an accretion disc would 
be a flare above the disc accompanied by a collapse down- 
wards (Fig. |ic|a) . The velocities we observe in V404 Cyg are 
larger than in solar flares, but these likely will rescale in a 
very different environment. There is a qualitative difference, 
however, in that our asymmetry persists well after the flare 
peak, and appears to be more pronounced on the decline 
than on the rise. 

Another mechanism emerges from the model of Haswell, 
Tajima & Sakai (1992). This model predicts explosive recon- 
nection of wound-up magnetic fields, accompanied by a bulk 
motion of plasma towards the reconnection point. In this 
case material closest to the observer will be red shifted and 
emission is expected over a range of velocities (Fig. |Io|b). 
The temporal properties might be expected to be similar to 
solar flares, however, with the asymmetry most prominent 
on the flare rise; this is not supported by our observations. 

Strong red peaks are actually not uncommon in 
BHXRTs, but are usually a feature of the outburst. Wu 
et al. (2001) have discussed the formation of lines in X- 
ray binaries and suggest that the asymmetric profiles seen 
in outburst are a consequence of absorption by a relatively 
slow disc wind driven by external irradiation. These asym- 
metric profiles are often accompanied by broad absorption 
troughs which are attributed to absorption by higher veloc- 
ity material. In one case, GRO J1655-40, we may also have 
seen P Cygni profiles in the ultraviolet, also indicative of a 
disc wind (Hynes et al. 1998) . In the case of V404 Cyg, we 
only see the asymmetric profile without the broad absorp- 
tion component. A plausible interpretation of this would be 
that a disc wind, or some other outflow, can be formed dur- 
ing strong flares. This could originate either in a flare in 
the central source which photoionises the disc, or in out- 
ward moving material from a reconnection event. In these 
flares, the whole profile is enhanced, but the blue wing is 
absorbed by approaching material (Fig. [To|c,cI) . No higher 
velocity ejecta are present, so there is no broad absorption 
component. Given that V404 Cyg remains a variable radio 
source in quiescence, it is also possible that there may be 
connection between such ejecta and the radio variability, 
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Figure 10. Possible geometries giving rise to a redward asymme- 
try during flares, a) A flare in the disc corona collapsing down- 
wards onto the disc. The emitting material is rcdshiftcd. b) A 
flare within the disc is accompanied by motion towards the re- 
connection point. The material nearest the observer is redshiftcd. 
c) A flare is accompanied by an outflow. The portion of the out- 
flow between the observer and the flare is blucshiftcd so absorbs 
the blue wing of the flare profile. The outflow could originate from 
a reconnection event or be associated with a flare of the central 
source, d) The line emission from the flare is spread over the 
whole disc. An outflow, either from the central source or the disc, 
partially absorbs the blue wing. 



with continued ejection of matter driving the radio variabil- 
ity. 



6.6 The quasi-periodicity 

The final question to be considered is the origin of the 6 hr 
quasi-periodicity of the flaring (Pavlenko et al. 1996). This 



was clear in earlier data, but less conclusive in 1999, al- 
though it may still be present. If flares originate in local 
magnetic reconnection events it is very hard to envisage a 
mechanism which would make events in different parts of 
the disc coordinate on a 6 hr timescale. 

Esin et al. (1997) suggested that interactions between 
an ADAF and the outer thin disc might produce quasi- 
periodic variations with a characteristic timescale of a mul- 
tiple of the Keplerian period at the transition radius. If we 
identify 6 hr with the Keplerian period at the transition ra- 
dius, r tr , then for a 12 M black hole this corresponds to 
logrtr ~ 4.8 in units of the Schwarzschild radius. This is 
much larger than assumed in ADAF models and indeed the 
high velocity wings of Ha rule out such a large transition 
radius (logrtr < 4.4; Narayan et al. 1997). There remains 
the possibility that the QPO period is a multiple (2 or 3) 
of the Keplerian period, although why this should be the 
case is not clear. Alternatively, 6hr could be some instabil- 
ity timescale intrinsic to the inner flow. 



7 CONCLUSION 

We have found significant spectral variability in the BHXRT 
V404 Cyg in quiescence; most nights exhibit some variabil- 
ity. The variations are most dramatic in Ha, (and presum- 
ably other emission lines as well), with changes of nearly a 
factor of two in the line strength over a few hours. Similar 
behaviour is seen in the continuum and this is usually cor- 
related with line changes. The continuum variations show 
smaller amplitude, as expected if the dominant continuum 
source, the companion star, does not vary. The variability 
can be phenomenologically divided into two types. Low level 
flickering involves changes in the line flux, without dramatic 
changes in the profile shape. Large flares seen approximately 
once per night, and may be associated with the ~ 6 hr QPO, 
involve the development of pronounced asymmetry in the 
line profile with the red wing strongest. 

We have considered several interpretations of the vari- 
ability. Chromospheric activity is seen on other rapidly ro- 
tating late-type stars and this gives rise to Ha flares. The 
energy released by the flares in V404 Cyg, however, is two 
orders of magnitude higher than in such systems, and the 
line profile changes are hard to reconcile with this interpre- 
tation so we consider it unlikely. Variability could originate 
from variations in the mass transfer rate from the compan- 
ion star and be manifested by flickering at the stream impact 
point. The kinematics of the observed flares do not appear 
consistent with this either, however. 

Variability could originate locally within the disc due 
to magnetic reconnection events, and may involve photoion- 
isation by X-ray flares in the corona immediately above the 
disc. To explain the participation of the whole line profile 
in the events requires a large spread in velocity of material 
involved in the flare, but this would be expected from the 
explosive reconnection model of Haswell et al. (1992). In 
this interpretation the redward asymmetry could indicate 
either motion of the flaring material away from us or par- 
tial absorption by material moving towards us. It is unclear, 
however, why the flare peak velocities do not vary according 
to where in the disc the flare originates. 

Instead the observed flare kinematics suggest that: i) 
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the center of light of the flares does not move significantly 
over an orbital cycle and ii) the profile of the flare is 
quite similar to the disc profile, except for the stronger 
red wing. The most obvious interpretation of is that the 
whole disc, centred on the slowly moving compact object 
(Ki ^ lSkms" 1 ), participates in the flare. It is hard to see 
how a magnetic reconnection event could affect the whole 
disc in this way, so it is most likely that the flares are in- 
duced by variable photoionisation from the X-ray source. We 
have compared the luminosity of the flares with the expected 
photoionising flux and find that this interpretation is plausi- 
ble; all advective-type models (ADAFs, ADAFs with winds, 
CDAFs) could produce the observed Ha emission, although 
models which predict most of their emission from large radii 
in the flow can more efficiently irradiate the disc. In this 
interpretation the redward asymmetry of the flare profiles 
can best be explained by absorption of the blue wing, either 
by an irradiation-driven disc wind or an outflow from the 
central source. The correlation between the line and contin- 
uum which we observe is expected for ADAF models which 
predict that the optical accretion continuum is self-absorbed 
synchrotron from near the compact object. In CDAF mod- 
els, continuum emission is attributed to the outer disc, but 
it could then result from reprocessed irradiation and so be 
correlated with variations in photoionised lines. The mecha- 
nism responsible for the 6hr characteristic timescale in the 
flaring remains unidentified. 

If this interpretation is correct then Ha flares should 
be accompanied by correlated X-ray flares; this can readily 
be tested with simultaneous observations. If such a correla- 
tion is not seen then a local flare hypothesis becomes more 
likely. A possible test of the magnetic recombination flare 
hypothesis, would be to measure linear polarisation in the 
Ha profiles during the impulsive phase of the flares since 
Ha an d H/3 have been found to be linearly polarised in solar 
flares ( Firstova fc Boulatov 1996 ), with a degree of polari- 
sation ~20 percent. 
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